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INVISCID FLOW ABOUT AIRFOILS AT HIGH SUPERSONIC SPEEDS 
By A. J. Eggers, Jr., and Clarence A. Syvertson 


SUMMARY 


Steady flow about curved airfoils is investigated analytically at 
higb supersonic speeds. Assuming air bebaves as an ideal diatomic gas, 
it is found tbat small pressure disturbances emanating from the surface 
of an airfoil are almost completely absorbed in the leading-edge shock 
wave (or a shock wave emanating from any other location on the surface), 
provided the flow deflection angles are not too close to those corre- 
sponding to shock detachment. This result is found to be essentia ll y 
independent of Mach number. As a consequence, it is shown that within 
the limitations of the assumption of .ideal gas flow, the shock -expans ion 
method may be used with good accuracy to predict pressure distributions 
on curved airfoils at arbitrarily high Mach numbers. This observation 
is verified with the aid of the method of characteristics applied to a 
10-percent-thick biconvex: airfoil at 0° angle of attack. It is further 
shown that the shock-expansion method can be easily employed to con- 
struct the entire flow field about a curved airfoil, accounting for 
shock-wave curvature and resulting entropy gradients in the flow. 

An approximation to the shock-expansion method for thin airfoils 
at high Mach numbers is also investigated, and is found to yield pres- 
sure distributions in error by less than 10 percent at Mach numbers 
above 3 and flow deflection angles up to 25°. This slender -airfoil 
method is relatively simple in form and thus may prove useful for some 
engineering purposes. To this end, tables are presented to facilitate 
its use. 

Effects of caloric imperfections of air manifest in disturbed flow 
fields at high Mach numbers are investigated, particular attention being 
given to the reduction of the ratio of specific heats from 1.4 toward 
1.0. So long as this ratio does not decrease appreciably below 1.3* it 
is indicated that the shock-expansion method, generalized to include 
effects of these imperfections, should be substantially as accurate as 
for ideal gas flows. This point is checked by comparing pressure distri- 
butions predicted by the' generalized shock-expansion method and a gen- 
eralized method of characteristics. Both methods are employed in forms 
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applicable for local air temperatures up to about 5000° Rankine, corre- 
sponding, for slender airfoils, to Mach, numbers up to the order of the 
so-called escape Mach number. Caloric imperfections caused reductions 
in the pressure coefficients below those predicted for flows of an 
ideal gas. In turn, there is- a general reduction in force and moment 
coefficients up to 10 percent. 

Hie slender-airfoil method is modified to employ an average value 
of the ratio of specific heats for a particular flow field. This sim- 
plified method has essentially the same accuracy for Imperfect gas 
flows as its counterpart has for ideal gas flows. 

An approximate flew analysis is made at extremely high Mach numbers 
where it is indicated that the ratio of specific heats may approach 
close to 1. It is found that 1 the shock-expansion method may be in con- 
siderable error as disturbances incident on the leading-edge shock wave 
are no longer largely absorbed in the wave. In this case, however, the 
Busemann method for the limit of infinite free -stream Mach number and 
specific heat ratio of 1 appears to apply with reasonable accuracy. 


nraionjciTON 


Small-disturbance, potential-flow theories have been employed 
widely, and for the most part successfully, for predicting the pressures 
(and velocities) at the surface of an airfoil in steady motion at low 
supersonic speeds. Thus the linear theory of Ackeret (reference l) has 
proven particularly useful in studying the flow about relatively thin, 
sharp-nosed airfoils at small angles of attack, while the second-order 
theory of Busemann (reference 2) has found application when thicker 
airfoils at larger angles of attack were under consideration. At high 
free -stream Mach numbers the range of applicability of any potential 
theory is seriously limited, however, due to the production of strong 
shocks by even the relatively small flow deflections caused by thin 
airfoils. The assumption of potential flow is invalidated, of course, 
by the pronounced entropy rises occurring through these shocks. 

This limitation on potential theories was early recognized and led 
to the adoption (see reference 3) of what is now commonly called the 
shock-expansion method. The latter method derives its advantage over 
potential theories princip all y by accounting for the entropy rise 
through the oblique shock emanating from the leading edge of a sharp- 
nosed airfoil. Consequently, so long as the disturbed air behaves 
essentially lik e an ideal gas, and so long as entropy gradients normal 
to the streamlines (due to curvature of the surface) do not signifi- 
cantly influence flow at the surface, the shock-expansion theory should 
predict the pressures at the surface of an airfoil with good accuracy - 
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it is tacitly assumed, of course, that the flow velocity is everywhere 
supersonic, and that the Reynolds number of the flow is sufficiently 
large to minimize viscous effects on surface pressures. 

The departure of the behavior of air from that of an ideal gas at 
the temperatures encountered in flight at high supersonic speeds has 
been the subject of some investigation in the case of flows through 
oblique shock waves. In reference 4, the effects of thermal and caloric 
imperfections on the pressure rise across an oblique shock wave was 
investigated at sea-level Mach numbers of 10 and 20 and it was found 
that these effects decreased the rise by less than 5 percent for maximu m 
temperatures up to 3000° R (corresponding to flow deflection angles up 
to 24°) . This decrease was found to be due almost entirely to caloric 
imperfections or changes in vibrational heat capacities of the air pass- 
ing through the shock wave. The changes in temperature and density of 
the air passing through the wave were affected to a considerably greater 
extent. Subsequently, an investigation was carried out by Ivey and 
Cline up to Mach numbers as high as 100 (reference 5) using the results 
for normal shock waves obtained by Bethe and Teller considering effects 
of dissociation (reference 6) . As would be expected, the pressures 
were found to be affected to a somewhat greater extent at the higher 
Mach numbers. 

The extent to which flow in the region of the leading edge of an 
airfoil departs from the simple Prandtl-Meyer type has also been inves- 
tigated at high supersonic airspeeds. If the surface is curved, for 
example, to give an expanding flow downstream of the leading edge, 
expansion waves from the surface will interact with the nose shock wave, 
thereby curving it and yielding a nonisentropic flow field. This flow 
field may be characterized not only by disturbances emanating from the 
surface but also by disturbances reflecting to some extent from the 
shock wave back toward the surface. The manner in which these phenomena 
dictate shock-wave curvature and surface pressure gradient in ideal gas 
flows at the leading edge has been treated by Crocco (reference 7) and 
more recently by Schaefer (reference 8)„ M unk and Prim (reference 9)j 
and others. In the cases considered by Munk and Prim it was found that 
surface pressure gradients were less (in absolute value) than those 
obtained assuming Prandtl-Meyer flow at the higher Mach numbers (i.e., 
Mach numbers greater than about 3) although generally by no more than 
about 10 percent. Since curved airfoils are likely to be of fundamen- 
tal interest at high flight speeds (see, e.g., reference 10), these 
phenomena would appear to merit further investigation, particularly as 
regards their influence on the whole flow field. In addition it would 
appear desirable to consider effects of gaseous imperfections through- 
out the field. 
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Such an investigation has therefore heen undertaken in the present 
report using the method of characteristics to accurately obtain flow 
fields, and as a basis for obtaining the more approximate methods of 
analysis. Bie method is employed in a generalized form -which allows 
caloric imperfections as well as entropy gradients in the flow to he 
considered at temperatures up to the order of 5000° R - thermal imper- 
fections are neglected as being unimportant in atmospheric air flows 
(see reference 4) . A 10-percent-thick biconvex airfoil is treated at 
Mach numbers from 3.5 to infinity, and the results are compared with 
the predictions of the shock -expans ion method, including a simplified 
form of the method applicable to slender airfoils at high Mach numbers, 
and a generalized form of the method including effects of caloric 
imperfections . 


a 
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c z 
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SYMBOIS 


local speed of sound, feet per second 
chord, feet 

characteristic coordinates (C-l positively inclined and C 2 
negatively inclined with respect to the local velocity 
vector) 

section drag coefficient 
section lift coefficient 

section moment coefficient (moment taken about leading edge) 


pressure coefficient 



specific heat at constant pressure, foot-pounds per slug °R 
specific heat at constant volume, foot-pounds per slug °R 
Mach number (ratio of local velocity to local speed of sound) 
static pressure, pounds per square foot 
dynamic pressure, pounds per square foot 
gas constant, foot-pounds per slug °R 
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s,n 

T 

t 

V 

x,y 

a 

P 

7 

5 

0 

P 

a 

u 

0 

A,B, 

C,D. 

1 
N 
S 


rectangular coordinates (in streamline direction and nor ma] 
to streamline direction, respectively) 

temperature, °E 

time, seconds 

resultant velocity, feet per second 
rectangular coordinates 


angle of attack, radians uni ess otherwise specified 


Mach angle, arc sine radians 


ratio of specific heats 



(Average value of y is 7 a .) 


flow deflection angle, radians unless otherwise specified 
molecular vibrational energy constant, °R (5500° B for air) 
mass density, slugs per cubic foot 


shock-wave angle, radians 


ray angle for Prandtl-Meyer flow, radians 


Subscripts 


free -stream conditions 

conditions at different points in flow field 
ideal gas quantities 

conditions just downstream of shock wave 
conditions on streamline 

Superscript 


vector quantities 
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DEVELOPMENT OF METHODS OF ANALYSIS 
Method of Characteristics 


Two-dimensional rotational supersonic flows have "been treated "by 
numerous authors with the aid of the method of characteristics, and 
various adaptations of the method have been found which are especially 
stilted for studying particular types of such flows. In the case of 
steady flows in which atmospheric air does not behave as an ideal dia- 
tomic gas, a very familiar and simple form of the compatibility equa- 
tions may be employed. To illustrate, consider the Euler equation 


P 


dV 

dt 


- grad p 


(1) 


the continuity equation 


div (pV) = 0 ( 2 ) 

and the equation for the speed of sound (evaluated at constant entropy) * 



(3) 


Rewriting equations (l) and (2) in the form of partial differential 
equations and transforming the resulting expressions to the character- 
istic or Ci, C 2 coordinate system, there is obtained, upon combination 
with equation ( 3 ), the following relations for steady flow: 


ctg P + = 0 

pV 2 VdCi Sc 2 / V^Ci Sc 2 / 

and 

c ~fcg P f frp , Sp V f & & \ = Q 

py 2 \Sci Bc 2 / \Sci Sc 2 / 


(b) 

(5) 


A simple addition or subtraction of equations (4) and (5) then yields 
the compatibility equations (see, e.g., reference U) 


- — = — p V 2 tan p 
0C1 


Ss 

Sc x 


(6) 
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and 


pV s tan 0 

^C 2 c)C 2 


(7) 


Now, in reference 4 "both caloric and thermal imperfections of air were 
considered and it was found that the latter imperfections 1 have a neg- 
ligible effect on shock processes in atmospheric air. It may easily be 
shown’ that this conclusion also applies to expansion processes, and for 
this reason caloric imperfections , only, are considered in detail in 
the present paper. These imperfections become significant in air at 
temperatures greater than about 800° R and first manifest themselves as 
changes in the vibrational heat capacities with temperature. Thus, the 
specific heats, cp and c v , and their ratio, 7 , for the gas also change. 
The equation of state remains, however, 

p = pRT (8) 

and the specific heats are still related to the gas constant by the 
expression 


Cp - cv = R (9) 

Furthermore, it readily follows from the differential energy equation 
and these expressions that the speed of sound is given by the simple 
relation 


a 2 = 7 RT (10) 

Combining equations ( 8 ) and (10) and noting that sin 0 = a/v there is 
then obtained 


pV 2 = 


7p 


sin 2 0 


(H) 


Hence, on combining this equation with equations (6) and' ( 7 ), it is 
apparent that the familiar compatibility equations 


Sp _ -27p 65 

BCj_ sin 2 0 


( 12 ) 


1 Thermal imperfections usually appear in the form of intermole cular 

forces and molecular-size effects, and may be accounted for with 
additional terms in the equation of state. 
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and 

8p 2yp 68 

8C 2 sin 2 3 8C 2 

also hold for the more general type of flow under consideration. These 
equations are basic, of course, to two-dimensional characteristics 
theory, and, as will he shown later, form a convenient starting point 
for developing simpler theories of two-dimensional supersonic flow. 

In order to apply equations (12) and (13) , it is evident that the 
manner in which y and 3 or M are connected to p or 8 must he known. 
Relations implicitly connecting these variables at temperatures up to 
the order of 5000° R nay he rea dily obtained from the results of refer- 
ence 4 by simply e limina t ing the terms therein accounting for thermal 
Inperfections. Thus we have as a function of the local static tempera- 
ture and free-stream conditions 


• (13) 



For isentropic flow along a streamline, the pressure is related to the 
temperature by the expression 


where 


p_ = * W 

PN * (T) 

7 ± 



(16) 



*(T) = 


(IT) 



2H 
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If there is a shock •wave in the flow, 2 in particular a nose or leading- 
edge shock, then the following additional relations obtained with 
equations (8), (10), and (15) and the conditions for continuity of flow 
and conservation of momentum along a stre amlin e through the shock are' 
also required: 


£o _ 1 
% 2 


T r 


(l+7|j!% 2 )— (l+7 0 M 0 2 )+ 

% 


J 


Tr 


n2 


(l+7l#% 2 ) _ ~ (1+7 oMo 2 ) 


+ b 


T r 


L R 


L R 


( 18 ) 


and 


sin 2 cr 


7 r %% 2 _ ± 

7 0 m 2 



tan. 8 


1 1 
*** 7qMq 2 

ORn/Po)-! 


(19) 


(20) 


Using the local static temperature as a parameter, the term 
2yp/sin 2 0 in equations (12) and ( 13 ) may now be evaluated with equa- 
tions (l^) through (17).' Equations (l8) through (20) define the initial 
conditions downstream of a leading-edge or other shock wave in the flow 
field. Thus, equations (12) through (20) provide all the information 
necessary to calculate the flow about an airfoil by means of the method : 
of characteristics. As described in detail in appendix A, the calcula- 
tion is of three general types; namely, (l) calculation of conditions at 
a point in the flow field between the shock and the surface; (2) calcula- 
tion of conditions at a point on the surface; and (3) calculation of 
conditions at a point just downstream of the shock. Case (l) entails 
the use of both compatibility equations, while case (2) entails the 
use of the compatibility equation for a second-f amily characteristic 


2 If there are no shock waves, then the subscript R in equation (16) 
can, of course, be replaced with the subscript o. 
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line in combination "with, the equation of the airfoil surface, and 
case (3) involves the compatibility equation for a first-family line 
in combination with the oblique shock equations. With the aid of the 
three general types of calculations the entire flow field about an 
airfoil can be built up numerically using a computing procedure working 
from the leading edge downstream. In cases where changes in the vibra- 
tional heat capacities with temperature are neglected, the calculations 
are of course simp li fied since y of the gas can be considered con- 
stant, and temperature, pressure, and density ratios are simply the 
ideal gas functions of Mach number. 


Shock -Expans ion Method 


General .- This method of calculating supersonic flow of an ideal 
gas at the surface of an airfoil is well known, entailing simply the 
calculation of flow at the nose with the oblique shock equations and 
flow downstream of the nose with the Prandtl-Meyer equations. Deter- 
mination of airfoil characteristics in this manner requires a small 
amount of time, of course, compared to that involved when the method of 
characteristics is used, hence the advantage of the former method. The 
questions arise, however, as to exactly what the simplifying assumptions 
■underlying the shock-expansion method are, and what form the method 
takes (for calculative purposes) when the gas displays varying vibra- 
tional heat capacities. 

The matter of simplifying assumptions may perhaps best be con- 
sidered by employing equations (12) and ( 13 ) , the basic c ompatfb ility 
equations. If these expressions are resolved into the streamline 
direction and combined, noting that 


and 


l ( Je. + J*e. n ) 

3s 2 cos p \ 3c x 3c 2 / 


35 

3s 


2 COS P \ 3C! 



there is then obtained the relation 


1 - 


1 + 


35 / 30 ^ 

ds/dCg \ 2 7p 36 
3s/3Cx / sin 2 p 3s 


( 21 ) 


( 22 ) 


3p 

3s 


36/3c 2 


( 23 ) 
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defining the gradient of p along s. If flow along streamlines down- 
stream of the nose is of the simple Prandtl -Meyer type, however, we 
have 


dp _ 2yp 85 
3s sin 2 p 6s 


(24) 


Hence it is evident that the requirement for this type of flow is 


ds/dCi 

ds/dC 2 


< < 1 


(25) 


Equation (25) is, of course, simply an approximate statement of a well- 
known property of Prandtl-Meyer flows; namely, that flow inclination 
angles are essentially constant along first-family Mach lines. It fol- 
lows from equation (12) that if equation (25) holds, then the pressures 
will also he essentially constant along these lines. It does not fol- 
low, however, that the Mach number will he constant, or for that matter 
that the first-family characteristic lines will he straight (as is the 
case for isentropic expansion flaws about a corner) . In fact, it may 
easily he shown that the Mach number gradient along Ci is proportional 
to the local entropy gradient normal to the streamlines, and that the 
Ci lines are curved according to the change in M. Thus we see that 
there is really only one basic assumption underlying the shock-expansion 
method; namely, disturbances incident on the nose shock (or for that 
matter any other shock) are consumed almost entirely in changing the 
direction of the shock. 3 Within the limitations of this assumption it 
is evident that the method provides a relatively simple means for cal- 
culating the whole flow field about an airfoil, including effects of 
shock-wave curvature (see appendix B) . In general, of course, the 
validity of this assumption can only he checked by comparison of calcu- 
lations using this method with those using the method of characteristics. 

The shock -expansion method for a calorically imperfect diatomic gas 
is readily deduced from the equations previously obtained. For example, 
flow conditions at the leading edge of an airfoil cam be evaluated with 
the oblique shock-wave expressions (equations (l8) through (20)) and the 
expression for conservation of energy (equation (15)). The variation 
of flow inclination angle with pressure along the surface is then 


3 It is interesting to note that the assumption of Thomas (reference 12) 
that pressure is a function only of flow deflection angle and entropy 
is equivalent to this assumption. It follows, of course, that the 
most general solution obtainable with Thomas’s series representation 
of the pressure is that given by' the shock -expans ion method. 
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obtained by graphical! y integrating equation (24); namely. 


Ss - S H = eiz££- dp ( 26 ) 

p s 

where the variables 7 , p, and f 3 are evaluated uaing equations (l4) 
through (17), employing the static temperature as a parameter. When 
extreme accuracy is not essential, this rather tedious calculation can 
be avoided and a relatively simple algebraic solution of the flow down- 
stream of the nose can be employed . 4 The details of this solution are 
presented in appendix C. In the special case of flow at high supersonic 
speeds about slender airfoils, the whole calculation becomes particularly 
simple and warrants special attention. 

Slender airfoils .- If it is assumed that the local surface slopes 
are small compared to 1 and in addition that the free -stream Mach number 
is large compared to 1, it follows that a and (3 are everywhere small 
compared to 1. In this case equation (24) takes on the approximate form 


dp 

dS 



ypM 


(27) 


Furthermore, if it is assumed that y is constant at an average value 
7 a for a particular flow field (this assumption appears reasonable 
since in the temperature range up to 5000° R the change in 7 is less 
than 10 percent as shown in reference 4), then the Mach number and 
pressure may be related by the simple expression 

7 -1 

M = %(il) 27a ( 28 ) 


Equations (27) and (28) combine to yield the differential equation 


-A- ( 2 -) 

4% \Pn/ 


(7 a +i) 



(29) 


which readily integrates (between IT and S) to the form 


4 The tabulated results of Noyes (reference 13) may also prove useful in 
this case for Mach numbers up to 3. 
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PS 


P 


N 




2/ a 

r a -i 


New denoting 


=f(M Q 6 N ) 


(30) 


(31) 


and 


=g(Mo5 N ) 

■^o 


(32) 


there is obtained from the oblique shock equations, simplified to con- 
form with this analysis. 


f(Mo5 N ) = 


Mq 2 ^n 2 - 1 


and 




g(Mo8 n ) = ~ (/a-!) 

7 a + 1 


(33) 


(3*) 


where 


7 +1 

MqON = Mq 5 n 


J\ + (zs^ 1Stfis y 


(35) 


With equations (30) through (35) the pressures on the surface of an 

airfoil may easily be obtained. In terms of pressure coefficient we 
haye 


°P = 


7aMo‘ 


fPNyps\ _ 1 _ 

V p o/\PN/ 


(36) 


or 


Cp = 

* 7aMo 2 


s(M d 6jj) 1-f (Mo8 N ) ^1 - ^ 


f>a_ 

7a -1 -1 


(37) 
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The advantage of these slender-airfoil expressions lies, of course, in 
their relative simplicity and thus the ease of calculation ■which is 
inherent to them. It may he noted in this regard that the functions 
f (MqSjj) and g(M 0 5jj) can he calculated once and for all with equations 
(33) , (34), an d (35), provided the variation of 7 a with Mq&^ is known. 

This calculation has heen carried out for a constant value of y equal 
to 1.4, and average values of y assuming T 0 = 500° R . 5 The results 
are presented in table I. 

It should also he noted that the s lender-airfoil expressions of 
the shock-expansion method satisfy the hypersonic similarity law for 
airfoils first deduced hy Tsien (reference l4 ) , 6 A necessary condition 
for the validity of "these expressions is thus satisfied; however, the 
accuracy of the shock-expansion method, whether for slender airfoils or 
otherwise, remains to he investigated. Such an investigation is now 
undertaken with the aid x of the method of characteristics. 


INVESTIGATION OF FLOW ABOUT AIRFOILS 
AND DISCUSSION OF RESULTS 


This study is divided into two parts: first, a consideration of the 
effects of Mach number assuming air behaves as an ideal diatomic gas; 
and second, a consideration of the combined effects of Mach number and 
gaseous inperfections, with principal emphasis in the latter regard 
placed on the caloric imperfections previously discussed. 


Ideal Gas Flows 


The effects of Mach number of primary interest here are, of course, 
those which result from interaction between the leading-edge (or other) 
shock wave and small disturbances originating on the surface of an air- 
foil. Some insight into the nature and extent of these effects can he 
obtained in the region just downstream of the shock wave without regard 


^For a given value of To, Tu, to the accuracy of this analysis, is the 
ideal gas function of Mo^N* Thus, knowing Tjj, can he determined. 


The average value of 7 used is 7 a = 7 a ( M o 5 N^ 


Wi. 

2 


6 This fact was employed hy Linnell (reference 15) to obtain an expression 
for pressure coefficient equivalent to equation ( 37 ) for the case of 
constant 7 , and to obtain explicit solutions for the lift, drag, and 
pitching -moment coefficients of several airfoils at hypersonic speeds. 
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(see equation (25)) 


for the shape of the airfoil producing the shock. To this end it is 

6&/6C1 6p/dCi , / w 

convenient to consider the ratio = (see equation (25)) 

6s/6c 2 

■which may he termed "the disturbance strength ratio" since in the region 
under consideration it is a measure of the ratio of strengths of distur- 
bances reflected from the shock wave to disturbances incident on the 
wave. This ratio may be evaluated with the expression 


d8/dc x 

a 8 /dC 2 


27m 

.sin 2 

27'tiPg 
.sin 2 


_ | \ 
35 * "s' 

+ -I ) 

as 


si 11 (f% + “ cr) 

sin (% - 6 jj + ff) 


Which is easily obtained with the compatibility equations and the 
oblique shock-wave equations as shown in appendix D. This calculation 
has been carried out for Mach numbers from 3*5 to 00 (7 W = 1.4) and 
flow deflection angles approaching those corresponding to shock detach- 
ment (i.e., Mjj ~ l) and the results are presented in figure 1. It is 
evident that except near 1, the ratio is small (in absolute value) 

compared to 1 throughout the entire range considered - this observation 
also applies, of course, at lower supersonic Mach numbers. Thus it is 
indicated that almost all of an incident disturbance is generally 
absorbed in the shock wave, provided the air behaves like an ideal dia- 
tomic gas. 7 This result is substantially the same, of course, as that 
which is assumed in deriving the shock-expansion method of calculating 
flows about airfoils, and therefore yields some credence in the method 
for high Mach number as well as low Mach number applications. 

As an over-all check on the shock-expansion method, surface pres- 
sure distributions calculated thereby are compared in figure 2 with 
those obtained with the method of characteristics for a 10-percent-thick 
biconvex airfoil (a = 0°) operating at free -stream Mach numbers of 3»5, 
5, 7«5, 10, 15, and oo. Predictions of the slender-airfoil approx ima tion 
to the former method for high supersonic speeds are also shown. There 
is no apparent difference between the pressure distributions given by 
the method of characteristics and the shock -expansion method up to a 


7 This result is contrary to that obtained by Lighthill (reference l6) 
who reports that for hypersonic flows, a disturbance is reflected 
from a shock wave with opposite sign but essentially undiminished 
strength. Lighthill’s conclusion appears to be based on an incorrect 
evaluation of his results for the case of very high Mach numbers. 
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Mach number of 10. At Mach numbers of 10, and above, however, the 
latter method predicts pressures which are slightly low downstream of 
the nose, becoming progressively lower with increasing Mach number. 

This result would be deduced from figure 1 where it is observed that, 
at the Mach n umb ers under consideration, expansion waves incident on 
the nose shock are reflected back toward the surface as compression 
waves of relatively small but increasing strength with increasing Mach 
number. The effect of these waves does not become pronounced even at 
infinite Mach n umb er (see fig. 2(f)) and the shock -expans ion method is 
thus substantiated as being a reliable simplified method for predicting 
the flow about airfoils at high supersonic speeds, again, so long as the 
air behaves as an ideal diatomic gas. The further simplified slender- 
airfoil method also appears to be a good approximation over the entire 
range of Mach numbers, 8 although, as would be expected from the assump- 
tions made in its development, it is in somewhat greater error than the 
shock-expansion method at the lower Mach numbers . 

The relative accuracy at high Mach numbers of the slender -airfoil 
method, linear and second-order potential theories may be seen in 
figure 3* As might be expected, the slender -airfoil method is more 
accurate than linear theory at both Mq = 5 and 15 , and more accurate 
than second-order theory at Mo = 15. It is perhaps surprising to note, 
however, that at the lower Mach number of 5 the slender -airfoil method 
is also somewhat superior to the second-order theory. 

The pressure distributions of figures 2 and 3 have been employed to 
calculate the zero-lift drag of the biconvex airfoil, and the results of 
these calculations, along with additional predictions of linear and 
second-order theory, are shown In figure 4. Predictions of the shock- 
expansion method are, of course, in best agreement with those of the 
method of characteristics; while the slender -airfoil method, although 
slightly less accurate than the shock-expansion method, is apparently 
superior to both linear and second -order theory at Mach numbers above 3. 

The preceding findings verify that so long as the disturbance 
strength ratio is small compared to 1, the flow along streamlines is 
essentially of the Prandtl-Meyer type. If we choose, on the basis of 

these findings, a maximum absolute value for of 0.06 (note 

^ S5/SC! &/dC s 

the maximum value of — — for the cases presented in fig. 2 was 

ctd/^C 2 

8 The hybrid expression for pressure coefficient obtained by Ivey and 
Cline (reference 5) gives reasonably good results also, although not 
as accurate as the slender-airfoil method at the higher Mach numbers 
under consideration. 
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approximately 0.06 at Mq = »), the region in -which the shock -expansion 
method is applicable can readily he obtained from figure 1. The upper 
boundary line of this region is shown in figure 5 end it is evident 
that it lies only slightly below (about 1° in general) the line corre- 
sponding to shock detachment given approximately by the = 1.0 line. 
Almost the entire region of completely supersonic (ideal gas) flow is 
then covered by the method. 9 (See shaded area of fig. 5.) 

The question naturally arises concerning the corresponding range 
of applicability of the slender -airfoil method. This question may be 
answered in part by comparing separately the predictions of the method 
for oblique shock flows and expansion flows with those of the exact 
oblique shock equations and Prandtl-Meyer equations. Such a comparison 
is shown in figure 6 in terms of the percentage error in the pressure 
coefficients predicted by the slender -airfoil method. As would be 
expected, this method does not exhibit good accuracy over the wide 
range of applicability of the shock-expansion method; however, it is 
indicated that it should predict pressure coefficients with less than 
10-percent error down to Mach numbers as low as 3 for airfoils producing 
flow deflections tip to as high as 25°. 

As a further check on the utility of the slender -airfoil method, 
the pressure coefficients on the 10 -percent -thick biconvex airfoil have 
been calculated with this method and the shock-expansion method at a 
Mach number of 10 and angles of attack up to about 30 0 . 10 The results 
of this calculation are shown in figure 7 (see fig. 2(d) for a = 0°) 
where it is seen that the agreement is reasonably good even at the 
highest angle of attack. This fact is reflected in figure 8 showing 
the force and moment coefficients for the airfoil as a function of angle 
of attack. Little difference is observed in the force coefficients as 
calculated by the two methods, while the moment coefficients display 
more pronounced but nevertheless small differences at the higher angles 
of attack. 


9 If it is required as by Rand (reference 17) that the entire flow field 
be of the true Prandtl-Meyer type (i.e., that all flow properties be 
constant along first-family Mach lines and not just 6 and p), then 
the range of applicability of the shock -expans ion method would be 
appreciably smaller. However, it has been shown that this restriction 
is not necessary. 

10 These conditions are within the range of applicability of the shock- 
expansion method as defined in figure 4; hence the use of the method 
as a base of comparison seems justified. Since the shocks -expansion 
method is far less tedious to apply than the method of characteristics, 
it will be employed as such a base in subsequent calculations whenever 
the conditions being investigated have been determined to be within 
its range of applicability. 
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From these and previous considerations, the ranges of applicability 
of the shock -expansion and slender -airfoil methods for ideal supersonic 
gas flows are reasonably well established. It r em ai ns now to determine 
the manner and extent to which gaseous imperfections in the flow at 
higher supersonic speeds may alter these ranges, and the reasons 
therefor. 


Imperfect Gas Flows 


As a first step toward investigating the effects of gaseous imper- 
fections on the high Mach number flows under consideration, it is con- 
venient to extend our consideration of the disturbance strength ratio 
ds/dCi 

It is recalled that when air exhibits a constant value of y 

bb/b C 2 

equal to 1.4 (the value for an ideal diatomic gas), the disturbance 
strength ratio is small at arbitrarily large Mach numbers, provided the 
flow deflection angles are not too close to those for shock detachment. 
One of the most Important effects of gaseous imperfections is, however, 
to decrease y of the disturbed air below this value due to the excita- 
tion of additional degrees of freedom (e.g., vibrational) in the mole- 
cules at the high temperatures encountered at high Mach numbers . Indeed, 
at arbitrarily high Mach numbers it might be expected that y of the 
disturbed air would approach 1 , since the number of degrees of freedom 
may effectively become very large (see, e.g., references 3 and 6 ). In 
this case, however, the extent of the disturbance flow field is decreased 
to a layer at the surface of the body which is negligibly thin compared 
to that for the case of ideal gas flow. Thus it is apparent that signi- 
ficant changes in the flow about airfoils at high Mach numbers may 
result from decreases in 7 of the disturbed air; hence the effects of 
such decreases on the disturbance strength ratio would appear to warrant 
attention. 

A detailed analysis of these effects is impractical at the present 
time due to the limited range over which the variation of 7 with tem- 
perature is known. However, some knowledge of these effects can be 
gained by repeating the ideal gas calculations for constant values of 
7 jj between 1.4 and 1.0 . 11 Such calculations have been carried out at 

1;L Since the enthalpy is negligibly small compared to the mass kinetic 
energy of the undisturbed fluid at the high Mach numbers of interest 
and, hence, 7 of this fluid does not influence the flow, this approach 
corresponds to employing an average value of 7 for the disturbed 
fluid. Since only flows of dense air are considered here, heat- 
capacity-lag phenomena are neglected (see references 5 and 6 ). 
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<3s/d Ci 

infinite Mach number since in this case has its may 1 mum value 

ds/dC 2 

for a given 7^, and the results are presented in figure 9* It is seen 
that except near shock detachment, the disturbance strength ratio 
increases with decreasing 7^, approaching 1 as 7^ approaches 1. This 

increase is slow at first; for example, the value of ^ P /^P 2 is still 

ds/dC 2 

less than 0.1 at 7-^ = 1.3. It might therefore be expected that the 
shock -expans ion method would continue to apply with reasonable accuracy 
so long as 7 of the disturbed flow is not appreciably less than this 
value . This point has been checked with the methods developed pre- 
viously for analyzing the flow of a calorieally imperfect diatomic gas 
at local air temperatures up to about 5000° R (note 7 has a value 
only slightly less than 1.3 at this temperature) . In particular, the 
pressure distribution on the lower surface of the biconvex airfoil at 
Mo = 10, <x=19.9°j and T 0 = 500° R (Tfj^l+OOO 0 R at leading edge) has 
been calculated with both the method of characteristics and the shock- 
expansion method. 12 The results of these calculations are presented in 
figure 10 and it would appear that the conclusions drawn from figure 9 
pertaining to cases where 7 N is of the order of 1.3 or greater are 
substantiated. Pressures in the expansion flow about the upper surface 
are not influenced (due to the low temperatures) by caloric imperfections 
and hence are the same as shown in figure 7(b). 

Shown also in figure 10 is the pressure distribution obtained by 
the shock -expans ion method for an ideal gas (7^ = 1.4). It is apparent, 
on comparing this pressure distribution with the other distributions, 
that although the effect of caloric imperfections on the disturbance 
strength ratio is small, the pressures are appreciably reduced by the 
increase in specific heats. The extent of this reduction is more com- 
pletely illustrated in figure 11 where the lower ' surf ace pressure distri- 
butions on the biconvex airfoil are presented for Mo=10 and T o =500 o R, 
at a = 0°, 10°, 19.9°> and 30°. As one might expect, the reduction in 
pressures increases with angle of attack (due to the corresponding 
increase in static temperature of the disturbed air) . The pressure 
coefficients calculated with consideration for the imperfections in the 
gas are less on the lower surface (up to 6 percent at the leading edge 
and 15 percent at the trailing edge) than those calculated assuming the 
gas behaves ideally. Hie upper-surface pressures are again unaffected 


12 For added ease of calculation the expansion method of appendix C was 
employed. This method is also employed in a l l subsequent calculations 
of this type since it has been found to yield results differing by 
less than 1 percent from those obtained by the more tedious graphical 
integration method. 
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by the caloric imperfections of air in all the cases presented (except 
at a = 0 °) since this surface experiences lower pressures and hence 
lower temperatures. They are therefore the same as shown in figure 7» 
Shown also in figure 11 are the pressure distributions calculated with 
the slender-airfoil method for y = y & . The accuracy of this simpli- 
fied method is substantially the same as was previously observed for 
the corresponding method in the case of ideal gas flows, although the 
local error may be greater than the reduction in pressure coefficients 
due to the caloric imperfections of air. This error is somewhat com- 
pensating, however, in its effects on the force and moment coefficients, 
as will be seen. 

The force and moment coefficients, corresponding to the lower- 
surface pressure distributions shown in figure 11 and the upper -surface 
distributions of figure 7 are presented in figure 12. The reduction in 
the lower-surface pressures leads, of course, to a general reduction in 
all three coefficients (up to about 10 percent for a = 30°) . The 
slender-airfoil method again predicts these coefficients with surprising 
accuracy. 

In order to further assess the accuracy, of the slender-airfoil 
method some additional calculations were carried out for the biconvex 
airfoil at a = 0° and Mq = 20 and 30. The pressure distributions for 
these cases were calculated by the shock-expansion method, slender- 
airfoil method (7 = 7 a ), and slender-airfoil method (7 = 7±) . These 
results are presented in figure 13 and it is observed that the use of 
7 a rather than 7 ^ Improves the accuracy of the slender-airfoil method. 
Tne extent of this improvement in the case of drag coefficient is shown 
in figure l4 - it would appear that predictions of the slender-airfoil 
method (7 = 7 a ) and shock-expansion method are in as good agreement as 
for ideal gas flows (see fig. 4) . On the basis of these and previous 
results, it may be concluded then that not only does the shock -expansion 
method retain its range of applicability when air exhibits caloric 
imperfections provided 7 of the disturbed air is not appreciably less 
than 1 . 3 , but also the slender-airfoil method (7 = 7 a ) retains its 
range. 


It would be surprising indeed, however, if this conclusion con- 
tinued to apply as 7 of the disturbed fluid approached 1 since, as 

Sd/aC-, 

discussed previously, is not small compared to 1 in this case, 

B5/Bc 2 

but would appear, in fact, to approach 1. This matter may be investi- 
gated in the same manner as the effect of 7 ^ on the disturbance 
strength ratio was investigated, namely, by using the ideal gas relation- 
ships in combination with appropriate values of 7 . 
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The limiting case of infinite free -stream Mach, n umb er and 7 = 1.0 
(for the disturbed fluid, see footnote 11 ) has already been investigated 
by Busemann (reference 18) and more recently by Ivey, Klunker, and Bowen 
(reference 19) . In this case, as pointed out previously, the shock wave 
e m anating from the leading edge remains attached to the surface down- 
stream of the leading edge (this is easily verified with the oblique 
shock— wave equations) and the disturbance flow field is co nfine d to an 
infinitesimally thin layer adjacent to the surface. In addition, the 
Telocity along a streamline downstream of the shock is constant, as may 
easily be shown with the compatibility equations. Surface pressures 
therefore become a simple function of airfoil geometry, 

_ d 8 o 

Cp = 2 sin 2 Sg + 2 cos 6 g S. 


varying, to a first approximation, directly with the square of the com- 
ponent of free -stream velocity normal to the surface (i.e., the flow is 
approximately of the Newtonian corpuscular type) . With this theory 
then, and the method of characteristics, we can get an idea of both the 
extent to which extreme changes of 7 from 1.4 toward 1 will alter 
surface pressures, and the accuracy with which the shock -expans ion 
theory predicts the alterations. To this end, figure 15 is presented 
showing the pressure distributions about the biconvex airfoil at 
Mo = oa as calculated by the several methods for different values of 7 . 
It is observed that, whereas the shock -expansion method agrees very 
closely with the method of characteristics for 7 = 1.4, there is a 
large difference at 7 = 1.05- This, of course, is precisely what one 
would expect from the previous discussion of the disturbance strength 
ratio. On the other hand, if the two characteristic solutions and the 
Busemann method are considered in order of decreasing 7 , it is indi- 
cated that the characteristics solutions approach the Busemann theory 
as 7 approaches 1. For 7 = 1.0 and Mq = “ the shock -expansion 
method, in turn, predicts a discontinuous pressure distribution with a 
pressure coefficient equal to that of the Busemann theory at the leading 
edge but a pressure coefficient of zero at al 1 points downstream of the 
leading edge. Hence it may be concluded that when the free -stream Mach 
number approaches infinity and 7 approaches 1, the Busemann method 
rather than the shock -expans ion method for calculating the flow about 
airfoils should be employed. 


I 


sin Sgdx 


(39) 


CONCLUDING REMARKS 


The flow about curved airfoils was investigated analytically at 
high supersonic speeds first assuming air behaves as an ideal gas, and 
then assuming air behaves as a thermally, although not necessarily 
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calorically, perfect gas. (Caloric imperfections had previously "been 
observed to be of predominant im portance in free flight, at least for 
local air temperatures up to the order of 5000° R.) It was found that 
so long as air exhibits no im perfections (i.e., behaves as an ideal 
diatomic gas) small disturbances originating on the curved surface of 
an airfoil are almost completely absorbed in the shock wave emanating 
from the leading edge (or any other location on the surface), provided 
the flow deflection angles are not too close to those corresponding to 
shock detachment. This result is essentially independent of Mach number, 
and is consistent with the early calculations of Crocco concerning sur- 
face pressure gradient at the leading edge of an airfoil. It was con- 
cluded that in ideal gas flows the shock -expansion method for determining 
flow conditions at the surface of an airfoil would apply with good 
accuracy at arbitrarily high Mach numbers. This conclusion was verified 
by the excellent agreement found between pressure distributions on a 
10 -percent -thick biconvex airfoil at 0° angle of attack calculated with 
the shock -expans ion method and the method of characteristics. It was 
further shown that the former method can be easily employed to construct 
the entire flow field about an airfoil in a manner that would account 
for shock -wave curvature and entropy gradients resulting therefrom. 

A high Mach number approximation to the shock -expansion method for 
thin airfoils was also investigated, and was found to apply with good 
accuracy at Mach numbers above 3 and flow deflection angles up to 25°. 

The essential feature of this slender -airfoil method is, of course, its 
simplicity, and for that reason it may prove useful for some engineering 
purposes. 

Effects of caloric imperfections were first investigated qualita- 
tively considering the reduction in the ratio of specific heats from 
1.4 toward 1.0. It was found that as the ratio decreased, the extent 
to which disturbances reflected from a shock wave increased. In the 
limit as the ratio approached 1, the reflection was complete, and the 
shock wave became tangent to the surface of the airfoil. So long as the 
ratio did not decrease appreciably below 1.3> however, it was found that 
less than 10 percent of a disturbance was reflected; hence the simple 
shock -expans ion method might be expected to continue to apply. This 
matter was checked quantitatively for the biconvex airfoil with the aid 
of a generalized method of characteristics including effects of caloric 
imperfections (up to local air temperatures of the order of 5000° R, 
corresponding to a ratio of specific heats of about 1.3). It was found 
that the shock -expans ion method was substantially as accurate as for 
ideal gas flows, provided it was also generalized to include effects of 
these imperfections. The principal effect of the reduction in specific 
heat ratio was to reduce the pressure coefficients below their ideal 
gas values by as much as 15 percent. The reduction in force and moment 
coefficients was somewhat smaller, being about 10 percent. Similarly 
it was found that the slender -airfoil method, modified to employ an 
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average value of the ratio of specific heats for a particular flow 
field, exhibited essentially the same accuracy as the analogous method 
for ideal gas flows. Thus it is indicated that the generalized shock- 
expansion method and its slender-airfoil counterpart can he applied 
with good accuracy up to very high free -flight Mach numbers. If the 
flow deflection angles are less than about 12°, these Mach numbers are 
of the order of the so-called escape Mach number (i.e., Mach numbers 
as high as 30 to 35) . 

At even higher Mach numbers -where the ratio of specific heats is 
expected to decrease appreciably below 1.3, and in fact perhaps to 
approach 1, it was not possible to obtain an accurate check on the pre- 
viously discussed qualitative considerations. It was undertaken, how- 
ever, to compare the pressure distributions on the biconvex airfoil 
predicted by the Busemann method (for the li mi t of the ratio of specific 
heats approaching 1 and Mach number approaching infinity) with those 
predicted by the shock -expans ion method and the method of characteris- 
tics at infinite Mach number (employing a constant ratio of specific 
heats of 1.05). As was expected, the shock -expans ion method was in 
very poor agreement with the method of characteristics, whereas the' 
Busemann method was in relatively good agreement. It is therefore 
indicated that for extremely high Mach numbers (something in excess of 
the sea-level escape Mach number) the Busemann theory may apply. 


Ames Aeronautical Laboratory 

national Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Jan. 9, 1952 
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APPENDIX A 

METHOD OF CHARACTERISTICS FOR TWO-DIMENSIONAL FLOW 
OF A CALORICALLY IMPERFECT GAS 


In the application of the method of characteristics for a c&lori- 
cally imperfect diat omi c gas to the particular problem of analyzing the 
flow about curved two-dimensional airfoils, many of the calculations 
are identical to those encountered in the solution of any problem where 
characteristics theory is employed. Since the details of these calcu- 
lations are well known, and well reported (see, e.g., reference 11), 
they will not he repeated here. 


A lattice -point system with an initial -value , numerical computing 
procedure will he used. The form of the compatibility equations to he 
enployed was developed previously; 13 however, it is convenient for pur- 
poses of calculation to substitute the pressure ratio, p/q^ into these 
equations and to rewrite them as difference equations. Equations (12) 
and ( 13 ) are thus reduced to the following forms 


and 


where 


(p/lo) c - (p/lo) A = ->^a(Sc - 5 A ) 

(Al) 

(p/qo) c - (p/Go) b = ( 6 C " 5 b) 

(A2) 

_ 27(p/qo) 

(A3) 

sin 2p 


It is also convenient to employ several. reference curves . These curves 
can he divided into two groups. The general reference curves consist of 
7 and ^(T) as a function of temperature, T. Equations (l4) and (IT) are 

13 This form of the compatibility equations (in p and 8 coordinates) 
was also used in obtaining some of the characteristics solutions for 
ideal gas flows. The majority of these solutions were carried out, 
however, with the compatibility equations in P, 8, and entropy 
coordinates , since it was found that greater accuracy was us uall y 
obtained for a given net size. In general, the net size employed 
yielded pressures at from 30 to 35 surface points on an airfoil with 
a maximum error in the corresponding pressure coefficients equal to 
less than 1 percent of the pressure coefficient at the leading edge. 
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used to determine these curves. A second set of shock -wave reference 
curves consisting of p /q Q , o', and 5 as a function of temperature, T, 
are determined hy use of equations (18) through (20) - the values of T 0 
and Mo are presumed known. 

In the computations three types of points are encountered. These 
are (l) a point in the flow field "between the shock wave and the air- 
foil surface, (2) a point on the airfoil surface, and (3) a point just 
downstream of the shock wave. Each one of these, types of points 
requires a slightly different computing procedure and they will be con- 
sidered in order. 


Point in the Flow Field Between' the Shock Wave 
and the Airfoil' Surface 


Figure 16(a) shows a schematic diagram of the system of points to 
he considered in these calculations. Point C is the unknown point at__ 
the intersection of the first-family characteristic line passing through 
point A and the second -family characteristic line passing through 
point B. Six quantities are known at both points A and B, and the 
problem is to calculate these same quantities at point C. These quanti- 
ties are x, y, 8, p/Iq, T, and %. The first five quantities are of 
obvious significance. The sixth, %, is defined as the static tempera- 
ture, just downstream of the shock wave, on the streamline passing 
through the point C. 

The physical coordinates of the point C (x^, y c ) may be determined 
by standard procedures such as those given in reference 11. In order to 
determine the quantity 8g, it is necessary to solve equations (Al) and 
(A2) simultaneously; thus 


6 C - 


- X A 5 A * X B 5 B + fe/q q) a ~ (p/1q) b 


X A + 




Equation (Al) or (A2) is then used to obtain (p/q,-,^ • 

There remains only the problem of determining Tq and Tjj at 

C 

point C. The temperature Tjy is obviously constant along the stream- 

C 

line through C. This quantity may therefore be calculated in the same 
manner as the entropy is calculated in similar flow fields for ideal gas 
processes (see, e.g., reference 11) .. Furthermore, since the flow along 
streamlines downstream of the shock wave is isentropic, equation (l6) 
may be applied in the form 
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(pAto) c _ »(?%) 

(p/qo)_ t(T c ) 


(A5) 


The pressure. 




7 


is defined in a manner analogous to 


%c » and 


may thus he determined using the shock -wave reference curves and the 
known value of Tjj- . Similarly, ijr ( T^ c ) may he determined from the 

general reference curves. The only unknown in equation (A5) then is 
■^(Tq) which may now he calculated. Once i); (Tq) is deter mine d, Tq may 
he determined hy again using the general reference curves. All six 
quantities, Xq, yg, 6g, (p/q^^ , Tq, and Tjj have now heen determined. 


Point on the Airfoil Surface 


Figure 16(h) shows a schematic diagram of the points to he con- 
sidered in these calculations. The physical coordinates of point C, 

(xq, y c ) , are first calculated hy solving simultaneously the equation 
of the second-family Mach line passing through point B and the equation 
of the airfoil surface. When xq and yg have heen determined, 8g is 
readily obtained from the equation of the airfoil surface. Equation (A2) 

is then applied to determine (p/qo) . 

C 

Since the airfoil surface is a streamline, Tjjg is constant along the 

surface and may he evaluated at the leading edge. The temperature, Tq, 

may then he determined using equation (A5) and the previously described 
procedure. All six quantities, xq, yg, (p/qo) c , 6g, Tq, and % c , are 
thus determined. 

In the special case of the first point on the airfoil surface 
downstream of the leading edge, the pressure ratio is calculated using 
the procedure of reference 9» This procedure is easily shown to he 
applicable to ealorically imperfect gas flows providing the oblique 
shock -wave equations of the present paper are employed. 

Point on the Shock Wave 


Figure 16(c) shows a schematic diagram of the points to he con- 
sidered in these calculations. The physical coordinates of point C 
(xq, Jq) are first calculated hy solving simultaneously the equation 
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of the first-family Mach line passing through point A and the equation 
of the shock wave linearized at point D, the last known point on the 
wave. The variation of p/q 0 with 6 along the shock wave may he 
approximated by the relation 


&/*>>C - ( ^D ■ ^ 


In this equation 


C 

a(p/qo) 

as 


(8p “ S-n) 


N 


is the rate of change of p/q^ 


N 


(AS) 

with 5 


along the downstream, side of the shock wave evaluated at point D. 

ie shock -wave equations, it is 

graphically or numerically 
N 

from the shock -wave reference curves. Equations (AL) and (A6) are 
solved simultaneously for Sq , thus 


Because of the complicated nature of ti 

generally easier to evaluate 

dS 


X A 5 A + 


cL(pAq) 


8„ = 


dS 


N 


8 d + (p/qo) - (p/Qq) 


D 


\a + 


d(p/qo) 


dS 


(A7) 


H 


When Sc has been calculated, Tq , and in turn (p/q^,)^ , may be deter- 
mined from the shock -wave reference curves. Since point C in this case 
is just downstream of the shock wave, Tq and T N are identical. 

c 

The six quantities, Xq, y c , (v/Iq) , 5 C , T c , and T^ have now been 

c c 

determined* 
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APPENDIX B 

SHOCK-EXPANSION METHOD FOR CALCULATING THE FLOW 
FIELD ABOUT AN AIRFOIL 


An initial value procedure -which, is similar to, although markedly 
simpler than, that associated, -with the method, of characteristics may 
"be employed, to carry out thi§ calculation. 14 ' To illustrate, consider 
the sketch: 



i 


With the oblique shock -wave and. expansion equations, all fluid, proper- 
ties at points M, A, C,. and. so forth on the airfoil surface may "be cal- 
culated. in the usuq.1 manner.- . If the point A is chosen close to M, the 
first-family (C x ) Mach line connecting A to point B on the shock wave 
may he con&idered-'straight and inclined at an ang^Le to the free -stream 
direction equal to 3 a + ^A*f Similarly, the segment MB of the shock 
wave may he considered straight and inclined at the angle cfy to the 
stream direction. Thus the ( point B in the flow field/may easily he 
determined. The direction of BD (a segment of the streamline passing 
through B) is the same as the. tangent to the surface at - A, and the 
attitude of the segment BE of the shock wave is fixed hy this direction. 
The locations of points' D and 'E in the flow field are thus fixed once 
point C is chosen. 15 The construction of the remainder of the flow field 
follows in a similar fashion. Having determined the shapes of the 
streamlines, the fluid properties along these lines are, of course, 
determined in the same mann er as those along the surface. 

14 It is clear, of course, that an "average value" procedure could also 
he employed. Such a procedure would, in fact, he the more desirable 
in some cases, since a coarser net may he used. 

The point C should, of course, he chosen close to A. 



NACA TN 26^6 


29 


It is important to note that this method, is applicable to the 
determination of the flow not only in the region adjacent to the air- 
foil (whether the -surface be concave or convex) but also in the region 
downstream of the airfoil; hence it may, for example, prove useful in 
downwash studies and the like. 


r 
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APPENDIX C 

APPROXIMATE SOLUTION FOR FRANDTL-MEYER FLOW 
OF A CALORICALLY IMPERFECT GAS 


The following solution is obtained with, an analysis sim il ar to 
that used in Meyer ’ s original paper (reference 20). A schematic dia- 
gram of the subject flow field is shown in figure 17 ♦ It is evident 
that the change in. flow-inclination angle for Prandtl -Meyer flow can 
he written as follows 


Sjj - 6 = (Pjj - p) + (wjj - oj) 


(Cl) 


Since the flow is isentropic, a given value of the local pressure will 
determine the Mach angle, p. The problem then is to evaluate the 
angle, cu. To this end the velocity components tangential and normal 
to the first-family Mach lines may be expressed in the usual maimer in 
terms of a potential op, thus 


Sep 

u = — 
Sr 


(C2) 


_ 1 ST 
r Sw 


(03) 


It is clear, however, that these components are functions of to only; 
hence it is convenient to define a new velocity potential which is a 
function of to alone. Such a potential is 

®( u ) = ^ (C40 


The velocity components may then be written in terms of this new 
potential. 


u = $ 

(05) 

a = A) 

(06) 


The resultant velocity is given by the expression 

V 2 = S’ 2 + G> u 2 


(07) 
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Equation (15) for 'conservation of energy may be written in terms of the 
local temperature as follows: 


V 2 + 2 


(A) 


]RT + 2RT ( — 

7± -lS \ e 0 /T_i. 


A c 


(08) 


The constant. A, is evaluated at the conditions existing upstream of the 
expansion region; namely, 


A * - + ( ^) ET » + 


(09) 


Equations (C7) and (c8) are then combined to yield 


<J> 2 + = -2RT 


0/T.V a2 
\r,-i ^e/T, y 


.7i-r e D/ ■ L -l 

It was shown previously, however, that 

a 2 = 7RT 


(CIO) 


(Cll) 


Equations (c6), (CIO), and (Cll) may therefore be combined to obtain 
the following relationship: 


0 2 + v 2 


1 + 


\ 7i -l e 0/T -l / J 


= A 2 


(C12) 


or 


/ 7 i +1 


7i 


1 7i-l + 7 ± - 1 L 7 


From the imperfect gas relationship for 7 we have 





1 + (7 ± -l) 


1 + 


7j-l 

7± 



e 0/T 

(e 0 /T-l) 2 


e 0 /? 

( e 8/T -D 2 . 


(Oli) 
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Substituting this relation into equation (C13) there is then obtained 

+ — F ( - ')]=A 2 (C15) 

l 7 ± -i y± V t / J 


< 5 2 + V 


where 


§) = 


0/T 


(e 0/T -l) 


1 -- 1, (7 1 (ee/T.p 


e/T r 


i+ 


0/T 


(e e ^-l) J 


1 + 


(7i-l) (6/T) 2 e 0 / T 
r i (e 0/T -l) 2 


(Cl6) 


Wow 


7 T = a 2 

?i 0 7 ± R0 


(C1T) 


For every value of T/0 there is thus a particular value of a 2 / 7jR0 . 
The function F(0/T) is therefore uniquely determined for any value 
of a 2 since 7jR0 is of course a constant. With this point in mind, 
let 


F(0/T) = G(a 2 /7 ± R0) (Cl8) 

Figure 18 shows G(a 2 /7-jR0 ) plotted as a function of a 2 / 7^R0 . This 
curve is approximated with the following simple relation: 


Gtaa/yjR 0) = 0.38 


7iR0 


+ 0.71 - 


O.lfr- 

2 /7iR0 


for 0.18 < -~ 
7 ± R0 

and 


< 1.0 


(C19) 


a 2 


G(a 2 /7 ± R0) = 0 


(C20) 
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Equation (C19) is also plotted in figure 18 to show the accuracy of this 
approximation. Consider first the case when G is given by equa- 
tion (C19) which is written in the form 

G(a2/7.Re) = g a 2 + p + t| \ (C21) 


\Jhe re obviously 


i = O.38/7 ± B0 

n = 0.71 

t] = -o.l 4 (7^0) 


(C22) 


Equation (C2l) is substituted into equation (C15) and with equations (Cl8) 
and (C6) the following expression results: 


— <t> 
7± u 


4 + 

Vri-i 




2 $2 4- ^ 3 . — * 2 


?i 


A" = 0 


(C23) 


In order to simplify this equation the following substitutions are made: 


and 


D 2 


( 


7j+l 

7i-l 


2eY + ± £ (A 2 — — T)) 

7±J 7 ± • 7± 


sin^v 



sin 2 T 



(C24) 

(C25) 

(C26) 


= 


D 2 


COS T 


(Tco) ; 


8 

?i 6 


( 027 ) 
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Equation (C23) then- reduces to 


, . u I cos V , ,2 4 Sin 2 3 t 4 sin 2 v ^ 

COS 4 5 6 T (Ty ) 4 + COS 2 T ( Ty) 2 + — = 0 ( C28) 


D c 


D 


This -equation is solved for t m , thus 


1 to 


= /- (cos T - COS v ) 1 / 2 

V D COS T 


(C2 9) 


or 


/ T) COS T dT 

du= /§ 


(cos T - COS v ) 1 / 2 


(C30) 


This expression is readily integrated to obtain the following equation 
relating to to the local velocity: 


io - tojjp = >j D ^2 
where 

E elliptic integral of the second kind 

F elliptic integral of the first kind 

k sin — (modulus) 

2 

z sin -1 ( sj ~ n (amplitude) 

V sin v /2 / 

The procedure for calculating corresponding values of the pres- 
sure, p, and the deflection angle, 5, is straightforward with the aid 
of the preceding equations and may he summarized as follows: 

1. Calculate A 2 , equation (C9) 

2. Calculate D 2 , equation (C24) 

3. Calculate V, equation (C25) 

4. Assume a value of T, less than T N 

5. Calculate p, equations (l 6 ) and (17) 

6 . Calculate V 2 , equation (C 8 ) 

7* Calculate 7 , equation (Cl4) 


E(k,z) - E(k,z N ) 


F(k,z) - F(k, z N ) 


(C3l) 
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8. Calculate a 2 (or $> u 2 ), equation (Oil) 

9. Calculate M and, in turn, 0 ftom 7 and a 

10. Calculate u 2 (or J> 2 ), equation (C7) 

11. CalciHate r, equation (C26) 

12. Calculate to , equation (C3l) 

13. Calculate 6, equation (Cl) 

This procedure is followed so long as the quantity a 2 / 7 -jR 0 is 
greater than 0.l8. (This is equivalent to the temperature being greater 
than approximately 1000° R.) For values of a 2 /y i R0 less than or equal 
to 0.18 (or temperatures less than about 1000° R) G is set equal to zero 
(see equation (C20)). In this case equation (C15) reduces to the same 
form as for an ideal gas, and therefore the well-known ideal gas rela- 
tionships can be used. 
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APPENDIX D 

EVALUATION OF THE DISTORBANCE STRENGTH RATIO 


Consider "the element of an oblique shock wave shown in the accom- 
panying sketch 



It is evident that the ratio 


ds/dCi 

is a measure of the ratio of 


d5/dC 2 

strengths of reflected and incident disturbances, respectively, since 
all disturbances incident on the _ wave between points D and C must 
travel along first -family characteristics which cross C s between 
D and A, while all disturbances reflected from the shock wave in this 
region must travel along second-family characteristics which cross C x 
between A and C. This ratio, termed the disturbance strength ratio, 
may be evaluated locally just downstream of the shock wave in the fol- 
lowing manner. The points D and C are chosen sufficiently close 
together so that the difference in pressure between these points may 
be written 


6p Sp 

Pc ' ^ “ ST AC2 + SE ACl 


or 


dp 

Pc-% — 


( 5 c -5 d) 


(EL) 


N 


(D2) 
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dp 

where — is the rate of change of pressure with flow deflection 

angle on the downstream side of the wave. Now the change in deflection 
angle between D and C is given by the relation 


R R 65 66 

= ^2 + 5cl AC 3 


Thus equations (D2) and (D3) combine to yield 


dp 

P C _P D = d5 


(l 


65 AH 

\ t ACo + 
n\^c 2 



(D3) 


(d4) 


But the compatibility equations (equations (12) and (13)) combine with 
equation (Dl) to give 


P C _I> D - 


sin 2Pjj 




(D5) 


Equating the right-hand members of this and the previous expression and 
rearranging, there is then obtained 


66/6C1 

65/6c 2 


( 2y W P N _ d£ N 

\ sin 2p N d6 w J 


’ 2y N% dp 
sin 2(3-^ d5 


) 



(D6) 


i^C2 

The ratio — — follows from the sine law, however: thus the disturbance 
AC X ’ ’ 

strength ratio is given by the relation 


65/6Ql 

65/6C 2 


( 2 ^ISPN dp 


\sin 2Pjj 36 

E / 


2y lffN dp 
in 2p N ^5 


N 


sin (Pfl+Sfl-g) 
sin (P|j~5jj+cr) 


(D7) 


which holds for both ideal and calorically. imperfect gas flows. 
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If the shock wave is attached to the airfoil at the leading edge 
(or more properly if >l), this expression may be used with equa- 
tion (23) to deter min e the surface pressure gradient at that point. A 
calculation analogous to this for ideal gas flows has already been car- 
ried out by Croeeo, Schaefer, Munk and Prim and others as discussed pre- 
viously. 
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TABLE I.- TABLE OF FUNCTIONS FOR SLENDER -AIRFOIL METHOD 


MoSjj- 

7 = 

= y± 

7 = 

7a (^0 - 

500°R) 



7 a 

g(Mo%) 


0 

1.000 

0 

1.400 

1.000 

0 

-05 

1.072 

.009901 

i.4oo 

1.072 

.009897 

.10 

1.148 

.01961 

i.4oo 

1.149 

•01959 

-15 

1.230 

.02912 

1.400 

1.230 

.02909 

.20 

1.316 

.03845 

l.4oo 

1.316 

.03841 

-25 

1.4o6 

.04760 

i.4oo 

i.4o6 

.04755 

-30 

1.502 

.05656 

1-399 

1.502 

.05649 

-35 

1.6o4 

.06537 

1-399 

1.603 

.06524 

.40 

1-710 

-07393 . 

1-399 

1.710 

.07380 

-45 

1.823 

.08235 

1-399 

1.822 

.08220 

-50 

1.94l 

.09058 

1-399 

i.94o 

.09040 

-55 

2.065 

.09863 

1-399 

2.064 

.09841 

.60 

2.195 

.1071 

1-399 

2.194 

.1062 

-65 

2.332 

.1142 

1-399 

2.330 

.1139 

-70 

2.474 

.1217 

1-399 

2.473 

.1213 

-75 

2.624 

.1290 

1-398 

2.622 

.1286 

.80 

2.780 

.1362 

1.398 

2.777 

.1356 

.85 

2.943 

.1431 

1.398 

2-939 

.1425 

-90 

3-112 

-14 99 

1.398 

3.108 

.1492 

-95 

3.289 

.1565 

1.398 

3.284 

.1557 

1-00 

3-473 

.1630 

1-397 

3-466 

.1620 

l.i 

3-862 

-1753 

1.396 

3.852 

.1740 

1.2 

4.280 

.1869 

1.396 

4.266 

.1853 

1-3 

4.728 

.1978 

1.395 

4.708 

.1958 

1.4 

5.206 

.2064 

1.393 

5.179 

.2056 

1-5 

5.715 

.2178 

1.392 

5.679 

.2147 

1.6 

6.256 

.2269 

1.391 

6.207 

.2232 

1-7 

6.827 

.2354 

1.389 

6.764 

.2309 

1.8 

7-431 

.2433 

1.388 

7.349 

.2381 

1-9 

8.066 

.2507 

1.386 

7.962 

.2447 

2.0 

8.734 

-2577 

1.384 

8.605 

.2508 

2.1 

9.411 

.2646 

1.382 

9.274 

.2564 

2.2 

10.17 

.2702 

1.380 

9.973 

.2616 

2.3 

10.93 

.2759 

1.378 

10.70 

.2664 

2.4 

11.73 

.2812 

1.376 

11.46 

.2709 

2.5 

12.56 

.2862 

1.374 

12.24 

.2749 

2.6 

13-42 

.2908 

1.372 

13.06 

.2788 

2.7 

14.32 

.2951 

1.370 

13.90 

.2823 

2.8 

15-25 

.2992 

1.369 

14.77 

.2856 

2.9 

16.21 

-3030 

1.367 

15.68 

.2887 
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TABLE I.- CONTINUED 


Mo% 

7 = 

7 ± 

7 = 7 a (To = 

500°R) 

g(Mo%) 

f(Mo8u) 

7 a 

g(Mo%) 

f(Mo%) 

3-0 

17-21 

0.3066 

1.365 

16.61 

0.2916 

3-2 

19.30 

-3130 

1.362 

18.56 

.2968 

3-4 

21.52 

-3187 

1.360 

20.64 

.3015 

3-6 

23.88 

.3237 

1.358 

22.83 

•3057 

3.8 

26.38 

.3280 

1.356 

25.16 

.3094 

4.0 

29.00 

.3322 

1.354 

27.60 

.3128 

4.2 

31.76 

.3358 

1.352 

30.17 

.3158 

4.4 

34.65 

.3389 

1.351 

32.86 

.3186 

4.6 

37.68 

-3418 

1.350 

35.68 

-3211 

4.8 

4o.84 

.3443 

1.349 

38.63 

.3234 

5-0 

44.14 

.3466 

1.349 

41.70 

.3255 

5.2 

47.56 

.3487 

1.348 

44.91 

.3275 

5-4 

51.13 

.3506 

1.347 

48.23 

.3292 

5.6 

54.82 

.3523 

1.347 

51.68 

.3309 

5-8 

58.66 

•3539 

1.346 

55-26 

-3323 

6.0 

62.62 

.3553 

1-346 

58.96 

• 3337 

6.2 

66.73 

-3566 

1-346 

62.79 

.3350 

6.4 

70.96 

-3578 

1.345 

66.77 

.3363 

6.6 

75-33 

.3589 




6.8 

79.83 

.3600 




7-0 

84.47 

.3610 




7.2 

89.24 

.3618 




7-4 

94.15 

.3626 




7-6 

99-19 

.3633 




7.8 

104.4 

.3640 




8.0 

109.7 

•3647 




8.5 

123.5 

.3661 




9-0 

138.2 

.3673 




9-5 

153-8 

.3684 




10.0 

170.2 

.3693 




10.5 

187.4 

.3701 




11.0 

205.4 

• 3707 




H.5 

224.3 

• 3714 




12.0 

244.1 

.3719 




12.5 

264.7 

• 3723 




13.0 

286.1 

.3727 




13.5 

308.3 

• 3731 




14.0 

331.4 

• 3735 




14.5 

355-4 

.3738 




15.0 

380.2 

.3740 
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TABLE I.- CONCLUDED 


Mq6 n 

•H 

i'- 

ll 

7 = 7 £ 

5. (T 0 = 500°R) 

g(Mo%) 



g(Mo&ff) 

f(Mo5 N ) 


16.0 

4-32.3 

0.37^5 




17.0 

487.7 

•37^9 




18.0 

546.5 

•3752 




19.0 

608.6 

• 3755 




20.0 

674.2 

• 3757 




22.0 

815.3 

.3761 




24.0 

969.8 

• 3764 




26.0 

1138 

-3766 




28.0 

1319 

• 3768 




30.0 

1514- 

•3770 




35-0 

2060 

• 3t72 




4o.o 

2690 

• 377^ 




45.0 

3404 

• 3775 




50.0 

4202 

• 3776 




60.0 

6050 

.3777 




70.0 

8233 

•3778 




80.0 

10750 

•3778 




90.0 

13610 

•3779 




100.0 

16800 

.3779 




00 

00 

.3780 
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Figure /,- Variation with deflection angle of the disturbance strength ratio behind an oblique shock 

wave for various free- stream Mach numbers (y e l.4). 
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(b) Mo =5 

Figure 2. - Pressure distribution on iO-percent-thick biconvex airfoil 
section for various free -stream Mach numbers at a =0.° 
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(d) Mo =10 

Figure 2. - continued. 



Pressure coefficient, C p Pressure coefficient, 
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Chordwise station, x/c 
(e) Mo = !5 



Chordwise station, x/c 
(f) Mo = oo 


Figure 2 \ — concluded. 
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Chordwise station, x/c 
(a) Mo =5 



Chordwise station, x/c 
(b) Mo = 15 


Figure 3.- Pressure distribution on /O- percent-thick biconvex airfoil 

section at a=0 ? 



Drag coefficient. 



Figure 4. — Variation of drag coefficient with free -stream Mach number for !0~ percent-thick biconvex 

airfoil section at a* 0° (y a i-4). 
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Free- stream Mach number, Mo 
(a) Oblique shock-wave flows 



Free- stream Mach number, Mo 
(b) Expansion flows 

Figure 6.— Accuracy of slender- airfoil method in predicting pressure 

coefficients (y=l4). 
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(b) a * 19.9 deg 



Chordwise station , x/c 
(c) a = 30.0 deg 



Figure 7. — Pressure distribution on 10 - percent- thick biconvex 
airfoil section for various angles of attack at Mo* !0 (7 s 14). 
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(a) Lift coefficient. 



(b) Drag coefficient. 



Angle of attack, a, deg 
(c) Moment coefficient. 


Figure 8.— Variation of force and moment coefficients with angle of attack 
for 10- percent- thick biconvex airfoil section at M» = 10 (y=l.4) 
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Calorically imperfect gas 
— — — Method of characteristics 

Shock-expansion method 

Ideal gas 

Shock -expansion method 
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Chordwise station, x/c 

Figure /O.- Pressure distribution on lower surface of !0~ percent-thick biconvex airfoil section at Mo* fO, 

To* 500° Ft, a = 19.9° 
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(a) a= 0° 




(b) a = 10.0° 


Figure II. — Pressure distribution on tower surface of iO -percent- thick biconvex 
airfoil section for various angles of attack at M, = iO and T, = 500° R. 
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Chordwise station, x/c 
(c) a=!9.9° 



Chordwise station, x/c 



(d) a = 30j0° 


Figure if. - Concluded. 
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Angle of attack, a, deg 

(a) Lift coefficient. 



Angie of attack, a, deg 

(b) Drag coefficient. 



Angie of attack, a , deg 


(c) Moment coefficient. 


Figure 12- Variation of force and moment coefficients with angle of attack 
for !0- percent- thick biconvex airfoil section at M, = !0, T, - 500° R 
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So 


MCA OT 2646 



Ghordwise station, 
(a) M.=20 



(b) M.=30 


Figure i3— Pressure distribution on !0- percent-thick biconvex airfoil section 


at (2=0 and Ti=500°R. 



Drag coefficient, c^ 


i 


t 



Free - stream Mach number, Mo 




Ftgure !4, - Variation of drag coefficient with 'Mach number for !0 ~ percent- thick biconvex airfoil section for 


H 


a=0° and To* 500° R. 
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Chord wise station, x/c 
(a) 1.4 



(b) x s i-05 


Figure t5. - Pressure distribution on iO~ percent- thick biconvex airfoil section 

at Mo -co and a - 0° 






9+I9S Mi VOVH 








